Reliable source of conditional non-Gaussian states from single-mode thermal fields 
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We address both theoretically and experimentally the generation of pulsed non-Gaussian states from classical 
Gaussian ones by means of conditional measurements. The setup relies on a beam splitter and a pair of linear 
photodetectors able to resolve up to tens of photons in the two outputs. We show the reliability of the setup 
and the good agreement with the theory for a single-mode thermal field entering the beam splitter and present a 
thorough characterization of the photon statistics of the conditional states. 



I. INTRODUCTION 

The subtraction of photons from an optical field is of both fundamental and practical interest, because it is linked to the 
properties of the annihilation operator and plays a leading role in quantum information protocols involving non-Gaussian states 
generation, manipulation and distillation. In fact, the simplest way to generate a non-Gaussian optical state starting from a 
Gaussian one consists in subtracting photons from it. Photon subtraction can be implemented by inserting a beam splitter in the 
optical path of the original state, detecting the number of photons of the reflected portion and selecting the transmitted portion 
only if a certain condition on the number of detected photons is satisfied. The challenging part of this scheme is the use of 
detectors able to resolve the number of photons. As a matter of fact, if, on one hand, it is nowadays quite easy to detect a single 
photon (see, e.g., Ref. [ 1] and references therein), on the other hand, the limited availability of photon counters that can resolve 
higher numbers of photons has led to the quest for indirect ways to obtain such information 

It is worth mentioning that the subtraction of photons allows not only the generation of non-Gaussian states, but also the 
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enhancement of the non-locality of bipartite states [5-9], or the generation of highly non-classical states II 1 OL II 1H useful for 



quantum informatio n purpo ses 11211 . Nevertheless, non-Gaussianity is a necessary ingredient for continuous-variable entan 



glement distillation [13 



1511 and different protocols relying on Gaussification of entangled non-Gaussian states Ilq - ll8ll or on 



de-Gaussification of entangled Gaussian states have been proposed 01911 . In all these approaches an important role is played by 
photodetectors able to perform conditional measurements. 

In this paper we report a thorough analysis of a setup based on hybrid photodetectors allowing the discrimination of the 



number of detected photons up to tens [ 20 



2 1H . The aim of the paper is twofold: firstly we demonstrate the feasibility of our 
setup and, secondly, we investigate its reliability by characterizing the generated conditional states. The input Gaussian states 
we employ to achieve these goals are single-mode thermal fields. Thermal states are diagonal in the photon-number basis, thus, 
the knowledge of their photon statistics fully characterizes them and their conditional non-Gaussian counterparts, which are still 
diagonal. Thanks to this property, we can give a complete analytical description of the behavior of our setup, including the actual 
expressions of the conditional states, and we can verify the agreement between the theoretical expectations and the experimental 
results with very high accuracy and control. This is a fundamental test in view of the application of our setup to the generation 



of more sophisticated states by conditioning non-classical, multipartite and multimode ones [22]. 

Throughout the paper we investigate two possible scenarios. We refer to the first one as "conclusive photon subtraction" 
(CPS): a photon-number resolving detector is used to condition the signal and to conclude which is the effective number of 
subtracted photons. The second one is the "inconclusive photon subtraction" (IPS): an "on/off" Geiger-like detector, i.e., a 
detector only able to distinguish the presence from the absence of photons is employed, preventing us from inferring the actual 
number of subtracted photons. 

The paper is structured as follows. Section [TT] addresses the generation of conditional states by means of detectors with an 
effective photon-number resolving power. We discuss the model in the presence of non-unit quantum efficiency and give some 
analytical results. In Section [III] we briefly review the IPS process on thermal Gaussian fields; we also investigate the main 
properties of the generated conditional non-Gaussian states that will turn to be useful for the characterization of our setup. In 
Section|lV]we report the experimental demonstration of our scheme and thoroughly characterize the obtained conditional states. 
Section[V]closes the paper and draws some concluding remarks. 



II. CONDITIONAL NON-GAUSSIAN STATES FROM THERMAL FIELDS VIA CONDITIONAL MEASUREMENTS 



In Fig.Q]we depict the conditional photon-subtraction scheme based on a beam splitter (BS) and two photon-number resolving 
detectors. Although in our experimental realization we will consider only thermal states, for the sake of generality here we 
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FIG. 1: Scheme for the generation of conditional non-Gaussian states via photon subtraction. A thermal input state g is mixed with the vacuum 
state go — 0) (0 at a beam splitter (BS) with transmissivity r. Two photon counters (R, T) with quantum efficiency rjk, k = R, T, are used to 
generate and analyze conditional states. See text for details. 

consider a diagonal state of the form g = ^ n g n \n)(n\. After the evolution through the BS with transmissivity t, the initial 
two-mode state Rq = g ® |0) (0| is transformed into the state 



R = E *n E Al(r)A-{r)\n -k)(n-l\® \k)(l\, 

n=0 k,l=0 



(1) 



where ^4" (r) = J (") r™ s ( 1 — r ) s . Then, the reflected part of the beam undergoes measurement. The positive-operator valued 



measure (POVM) describing a realistic photon counting device with quantum efficiency 77 is given by [23] 
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in which B s , m {rf) = (m) 7 ?" 7 ^! — r l) S m ■ ^ ^ e photon-counter in the reflected beam detects mpj photons, the conditional photon 
subtracted (CPS) state obtained in the transmitted beam is: 
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where r]k is the quantum efficiency of the detector located in the reflected (fc = R) and in the transmitted (fc = T) beam paths, 
respectively. Note that the state in Eq. (0 is still diagonal. The overall probability pr(to r ) of measuring ran m the reflected 
beam reads: 
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(4) 



which represents the marginal distribution of the joint probability, 



Ftr(to t ,tor) = Tr[i?II mT (?7 T ) <S) U mR (r)n)] 
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= B s,m R (r]R) Bt, mT (rn) 



x [Al +t (r)] 2 g s+t , (5) 
that detectors T and R measure mx and photons, respectively. By taking g in a single-mode thermal state u(N t h), 

oo 

v(N th ) = J2M N th)\n)(n\, (6) 

n=0 

^th) = — -rH . (7) 



where jV t h denotes the mean number of thermal photons, Eq. (0 reduces to 



m R 

PTR(rn T ,m R )=["" ± '_"'"■) ■ x ■ 'T mn , , m „ , , , (8) 
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Mt — ttjt N t h and Mr = (1 — r) t^r iV t h being the mean numbers of detected photons of the transmitted and reflected beams, 
respectively. 

Given tor and pr(tor), the conditional state £>cps m Eq. ((3) can be obtained straightforwardly. From Eq. (0 we can then 
evaluate the Fano factor 

_ a 2 (M CPS ) 

FCPS " ^cp^' (9) 
which is the ratio between the variance er 2 (Afcps) an d the mean number Mcps of the photons detected in the CPS state. As we 
will see below, Ft > i*bps > 1, where Ft = 1 + Mr is the Fano factor of the single-mode thermal field of the (unconditional) 
transmitted beam. Note that qcps is always super-Poissonian, which is consistent with the classical nature of the field. 

To deeply characterize the output conditional state we evaluate its non-Gaussianity. Since the state has the form pcps = 



J2 n Pn\ n )( n \> the non-Gaussianity (nonG) measure 12511 can be written as: 

$[qcps] = S[u(N CP s)} +^p„logp„ (10) 

n 

where Nqps is the mean photon number of £>cps> and S[v(N)] = N log(l + l/N) + log(l + N) is the entropy of the thermal 
state v{N). 

However, due to the inefficient detection, we cannot reconstruct the actual photon number distribution p n , but only the detected 
photon number distribution q mT = j?tr(tot, tor) given in Eq. ((H), where tor is the conditioning value, and tot is the number 
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of detected photons. Thus we can evaluate the quantity 

elffcps] = SMMcps)] + ^2 9mT 1o S*™t < *[ecps]- (H) 

The last inequality follows from the fact that the inefficient detection may be described by a Gaussian lossy channel that does 
not increase the non-Gaussianity, followed by an ideal (i.e., unit quantum efficiency) detection (see Appendix lAl for details). 
The quantity e[g>cps]> which can be easily evaluated from our experimental data, turns out to be a lower bound for the actual 
non-Gaussianity, that is, significant values of e[gcps] correspond to more markedly non-Gaussian states Qcps- 



III. INCONCLUSIVE PHOTON SUBTRACTION ON THERMAL STATES 



The conditional states introduced in the previous Section can be generated only if the detector in the reflected beam path is 
able to resolve the number of incoming photons. In this Section we consider a scenario in which the detector R (see Fig.Q]) can 
only distinguish the presence from the absence of light (Geiger-like detector): we will refer to this measurement as inconclusive, 
as it does not resolve the number of the detected photons. When the detector clicks, an unknown number of photons is subtracted 
from g and we obtain the IPS state pips- To characterize this class of conditional state, we use the phase-space description of 
the system evolution, that allows a simpler analysis with respect to that based on the photon number basis. 

The phase-space description of the IPS operated on single-mode Gaussian states can be obtained by generalizing the analysis 



given in Ref. Ill ill . The Wigner function of the thermal state in Eq. © reads as follows (in Cartesian notation): 

expf-iX^-;^ 

Wth(X) = 



where: 
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(13) 



is the covariance matrix (CM), 1 being the 2x2 identity matrix. According to ill ill , the action of the BS transforms the CM of 
the two-mode input state (thermal+vacuum) 

/ 

(14) 



as follows [24]: 
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where A, B, and C are 2x2 matrices and 



(16) 



is the symplectic transformation associated with the evolution operator Ubs of the BS. 



The probability p on = p on {r, t, ,) that the on/off detector endowed with quantum efficiency 77R clicks is given by [24]: 



Pan = 1 -PoB(r,T,TfR.) 



1 - [ri R y/Det[B + Ml 
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(17) 
(18) 
(19) 



where p a g is the probability of a non-click event and 



0"M = -7: 1. 

2?7r 



The Wigner function associated with the IPS state pips reads: 

W a (X)-p off W b (X) 



(20) 



W IPS (X) = 



where 
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(21) 
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2 7 rv / Det[S^ 

S a = A and St = A — C(J3 — er M )C T . Note that the IPS, being it the linear combination of two Gaussian functions, is no 

r~| 

longer Gaussian: for this reason the IPS process is also referred to as a de-Gaussification process [ 10]. The Wigner functions in 
Eq. (1221 are those of two thermal states v{N k ) with mean number of photons given by 

rN th 



N a = rN th , N b = - 

1 + Vr(1 - rjJVth 

respectively; thus, the density matrix associated with (fJTJ can be written as: 

v(N a ) - PoS v{N b ) 



(23) 



QIPS 



Po 



(24) 



and the corresponding conditional distribution of the detected photons is: 

{M a ) - PoS Vm T {Mb) 



p T {m T ) 



Po 



(25) 



where M a — rj^Na and Mb — ri^Nb, being rj^ the quantum efficiency of the photon-resolving detector of the IPS state, and the 
f{Nk) are given by Eq. (Q). 
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Starting from the above results, we can give some further detail about the IPS thermal state in Eq. (l24l i. The mean number of 
detected photons is 

Mips = M^M, 

Pan 



and the variance ct 2 (Mips) is: 



2, w ^ M a(l + M a ) - PoS M b (l + M b ) 
a (.-Mips J = 
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(27) 



Moreover, as M a > M b , the Fano factor F\ps is: 

= ^(Mzps) 
Afips 

MJM a - M 6 ) M a - M 6 
M a - p oS M b 1 - p oS 

in which the final inequality can be checked by substituting the actual expressions of M a , M b and p g. The state is always 
super-Poissonian (also in this case, as one would expect, Ft > -Fips > !)■ As for the conditional states qcps, a ls° m this case 
we can characterize the non-Gaussianity of the state pips from the experimental data by evaluating the quantity 

e[erps] = S(v{Mips)) + ^ qm ^ lo § qmr ? > ^ 30 ^ 
n it 

which is still a lower bound for the non-Gaussianity measure, i.e. e[gips] < <5[f?rps]> as explained in the previous Section and in 
the AppendixlAlin more detail. 



IV. RELIABLE SOURCE OF NON-GAUSSIAN STATES 
A. Experimental setup 

We produced a single-mode pseudo-thermal state by inserting a rotating ground glass plate in the pathway of a coherent field, 
followed by a pin-hole to select a single coherence area in the far-field speckle pattern (see Fig. |2}. As most detectors have the 
maximum quantum efficiency in the visible spectral range, we chose to exploit the second-harmonic linearly polarized pulses 
(A = 523 nm, 5.4 ps pulse duration) of a Nd:YLF mode-locked laser amplified at 500 Hz. The thermal light was split into 
two parts by a polarizing cube beam-splitter (PBS) whose transmissivity r can be continuously varied by means of a half-wave 
plate (A/2 in Fig. [2j. We balanced the two exiting arms of the PBS to achieve r ~ 0.5. The light exiting the PBS was focused 
in two multi-mode fibers and delivered to two hybrid photodetectors (HPDr .t, mod. R10467U-40 Hamamatsu), endowed not 
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FIG. 2: (Color online) Scheme of the experimental setup: P, rotating ground glass plate; PH, pin-hole; ND, continuously variable neutral density 
filter; A/2, half-wave plate; PBS, polarizing cube beam-splitter; L, collective lens; MF, multi-mode fiber; HPDr,t, hybrid photodetector; AMP, 
preamplifier plus amplifier; SGI, synchronous-gated integrator; ADC+PC, analog-to-digital converter. 

only with a partial photon-resolving capability, but also with a linear response up to 100 incident photons. The outputs of the 
detectors were amplified (preamplifier A250 plus amplifier A275, Amptek), synchronously integrated (SGI, SR250, Stanford), 
digitized (ATMIO-16E-1, National Instruments) and processed offline. To analyze the outputs we model the detection process 
as a Bernoullian convolution and the overall amplification/conversion process through a very precise constant factor 7, which 



allows the shot-by-shot detector output to be converted into a number of detected photons 12711 . The calibration procedure 



required performing a set of measurements of the light at different values of the overall detection efficiency of the apparatus, 
77, set by rotating a continuously variable neutral density filter wheel placed in front of the A/2 plate. For each value of 77, 
we recorded the data from 30 000 subsequent laser shots. For the results presented in the following we obtained the values 
7r = 0.104 V and 7t = 0.093 V for the calibration of the detection chains in the reflected and transmitted arms of the beam 
splitter, respectively. These values of 7r,t were used to convert the voltages into number of detected photons that were finally 
re-binned into unitary bins to obtain probability distributions. Once checked the reliability of the calibrations from the quality 
of these distributions, the voltage outputs of the HPDr.t detectors were associated with numbers ihr and mx in real time. The 
linearity of the detectors and the absence of significant dark counts make our system suitable for making experiments in both the 
CPS and IPS scenarios to produce conditional states. In the case of IPS, we only distinguish the HPDr outputs giving tor = 
from those giving any tor > 1 to mimic the behavior of a Geiger-like detector. 



B. Conditional non-Gaussian states 

The good photon-resolving capability of HPD detectors and their linearity make it possible to implement the CPS scheme de- 
scribed in Section|IIl Conditional measurements in the reflected beam irreversibly modify the states measured in the transmitted 
arm and in particular make them non-Gaussian. 
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FIG. 3: (Color online) Joint probability ptr(tit, ttir) to measure ttir photons in the reflected beam and ttit photons in the transmitted one. 
The experimental data (red dots) are plotted together with the theoretical surface (gray mesh). The other involved experimental values are 
Mr = 1.679 and M T = 1.254. 



To better understand the power and the limits of this kind of conditioning operation, we follow two different approaches: first 
of all we fix the energy Nth of the initial thermal field and characterize the CPS state as a function of the conditioning value tor. 
Secondly, we consider the properties of CPS states as a function of the mean incoming photons Nth for a particular choice of the 
number tor of photons detected in the reflected arm. The final aim is the production of non-Gaussian states with well defined 
conditioning value. 

We start by presenting the results obtained by choosing a set of measurements having Mr ~ 1.254. 

The joint probability ptr("t-t, wr) of measuring tor photons in the reflected arm and tot photons in the transmitted one 
is plotted in Fig. [3] as dots together with the theoretical surface to which they perfectly superimpose. Of course, starting from 
the theoretical joint probability, we can calculate the expected photon-number distribution of the states obtained by performing 
different conditional measurements in the reflected arm [see Eq. (0] and thus evaluate all the quantities necessary to characterize 
the CPS states. 

In Fig. |4] we plot the behavior of the mean number of photons Mcps of the conditional states and their Fano factors -Fcps as 
a function of the different conditioning values tor. We find that the Fano factor does not depend on the particular choice of the 
conditioning value tor, in agreement with the analytical result calculated from Eq. ((H): 

1 + M T + Mr 



F, 



CPS 



1.468. 



1 + Mr 

Note that the obtained value is definitely lower than that of the unconditional state, Ft 2.225. 



(31) 
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FIG. 4: (Color online) Fano factors Fqps of the conditional states: experiment (black dots) and theory (solid black line) as functions of the 
conditioning value ttir. The green line corresponds to the Fano factor Ft = 1 + Mr of the (unconditional) transmitted state. The inset 
refers to mean number of detected photons Mcps of the CPS states as a function of the conditioning value ma: experimental data (dots) and 
theoretical curve (solid line). The values of the other involved parameters are Mr = 1.679 and Mt = 1.254. 




FIG. 5: (Color online) Reconstructed photon-number distributions of the (unconditional) thermal state with mean number of photons Mt = 
1.254 (black triangles) and of the conditional states for six different conditioning values mn (colored dots and squares). The theoretical curves 
are plotted as lines according to the same choice of colors. The corresponding fidelity / is also reported. 

The photon-number distributions of the conditional states look quite different from each other. As it is shown in Fig. [5] the 
larger the conditioning value, the more different the statistics of the conditional state (colored symbols + lines) is from that 
of the incoming one (black triangles + dashed line). We note that, due to the limited number of recorded shots (only 30 000 
laser-shots), the experimental points tend to deviate from the expected px distributions at increasing conditioning values. This 
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FIG. 6: Lower bound e[f>cps] for the nonG measure <5[£>cps] as a function of the conditioning value tub, for Mt = 1.254. 



behavior can be quantified by calculating the fidelity (see / values reported in Fig. [5]): / = Ylm=o \j Pr 1 ( m ) Pr(™)> in which 
(m) and Pr{m) are the theoretical and experimental distributions, respectively, and the sum is extended up to the maximum 
detected photon number, m, above which bothp^(m) and pt(iti) become negligible. For all data displayed in Fig.[5]the fidelity 
is rather satisfactory. 

Finally, the behavior of the lower bound for the nonG measure as a function of the conditioning value tor (Fig|6|l predicted 
by the theory (line) is well reproduced by the experimental data (dots). In particular, it is worth noting that the value of e[gcps] 
increases at increasing the conditioning value. 

As an example of the second approach, we consider the CPS states obtained by choosing wir = 2 as the conditioning value 
for different values of Nth- In the inset of Fig.|7]the mean number of photons of the CPS states is plotted together with the mean 
number of photons of the initial states measured in the transmitted arm: it is interesting to notice that the values of Mcps actually 
approach the conditioning value tor = 2. Again, the experimental results (dots) are well superimposed to the theoretical curves, 
calculated starting from Eq. (0 with the measured mean values. 

Figure|7]also shows the comparison between the Fano factor of the unconditional states Ft (green squares) and that of the CPS 
states i^cps (black dots): as expected from theory the conditional states preserve the super-Poissonian nature of the incoming 
states, though with a smaller value of the Fano factor (Ft > -Fops > !)■ 

In Fig. |8]we show three examples of conditional-state photon distributions for different values of the total incident intensity. 
For each example, we plot both the original thermal distribution (full symbols) and that of the conditional state (empty symbols). 
The agreement with the corresponding theoretical predictions (lines) is again testified by the high value of the fidelities. 

We plot in Fig. |9]the lower bound for the non-Gaussianity s [gcps] as a function of the total mean detected photons (see Fig. [9]) 
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FIG. 7: (Color online) Log-Linear plot of the Fano factors Fcps of the CPS states (black dots) and of the unconditional states Ft (green 
squares) as functions of the total mean detected photons A/t + A/r for conditioning mrt = 2. The solid lines refer to the corresponding 
theoretical curves. The inset shows the mean number of detected photons Mqps of the CPS states as a function of the mean number of 
detected photons A/t of the unconditional states: experimental data (black dots) and theoretical curve (solid line). The green line refers to the 
mean photon number A/t of the unconditional states. 

together with the expected theoretical results. 

C. IPS non-Gaussian states 

Here we consider the scenario in which an on/off Geiger-like detector measures the reflected part of the input signal. In 
particular, as described in Section[III] we are interested in studying the properties of the state produced in the transmitted arm of 
the PBS whenever the detector placed in the reflected arm clicks. To this aim, we performed a set of measurements by fixing the 
transmissivity of the PBS r = 0.5 and changing the mean intensity of the light impinging on the PBS. 

In the inset of Fig. [10] we plot the mean number of detected photons Mips of the IPS states as a function of the mean 
number of detected photons Mr of the unconditional thermal states (black dots) together with the theoretical prediction (solid 
line) according to Eq. d26i >. We note that the effect of the conditioning operation is to increase the mean value of the original 
state. As described in SectionUm another quantity to characterize the IPS state is the Fano factor Fipg: to better appreciate the 
difference between the unconditional states and the corresponding conditional ones, we plot in the same figure (see Fig.flOb the 
corresponding Fano factors as functions of the total mean detected photons (symbols). For each set of experimental results we 
also plot the theoretical behaviors (lines): analogously to the conditional case, we have Ft > Fips > 1. 
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FIG. 8: (Color online) Reconstructed photon-number distributions for three different examples (black, green, red) of unconditional states (full 
symbols) and for the corresponding CPS states (empty symbols) with mR = 2. The theoretical curves are plotted as lines according to the 
same choice of colors. The corresponding fidelity / is also reported. 
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FIG. 9: Log-linear plot of the lower bound e[gcps] for the non-Gaussianity <5[f?rps] as a function of the total mean detected photons Mt + Mr 
and for tur = 2: experimental data (dots) and theoretical curve (solid line). 

In Fig.[TT|we show the reconstruction of the detected photons distribution Pt(tot) of both the unconditional (full symbols) 
and the conditional states (empty symbols) for three different mean values (black, green, red) of the incident intensity. The 
agreement with the corresponding theoretical distributions (colored lines), calculated with the measured mean values, can be 
checked by evaluating the fidelity, as reported in Fig.fTTl 

Finally, in Fig.[T2]we plot the lower bound for the nonG measure e[gips] as a function of the total mean detected photons. The 
correspondence between the experimental results (dots) and the theoretical prediction (line) is good. Note that e[grps] increases 
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FIG. 10: (Color online) Log-linear plot of the Fano factors fips of the IPS states (black dots) and of the unconditional states (green squares) 
as functions of the total mean detected photons Air + Mr. The solid lines refer to the corresponding theoretical curves. The inset shows the 
mean number of detected photons Mips of the IPS states as a function of the mean number of detected photons Mt of the unconditional states: 
experimental data (black dots) and theoretical curve (solid line). The green line refers to the mean photon number Mr of the unconditional 
states. 

as the mean number of detected photons increases: this allows the generation of highly populated non-Gaussian states. 



V. CONCLUDING REMARKS 

In this paper we have discussed in detail, both from a theoretical and an experimental point of view, a setup based on a single 
beam splitter and two photon-number resolving detectors to subtract photons from an incoming state and, thus, to generate 
non-Gaussian states starting from Gaussian ones. In order to show the reliability of our setup, we used (Gaussian) thermal states 
as inputs and completely characterized the conditional non-Gaussian outgoing states. In our analysis, we have adopted two 
possible detection schemes: the first one is based on the conclusive photon subtraction (CPS), whereas the second one on the 
inconclusive photon subtraction. In particular, we have demonstrated, as one may expect, that the non-Gaussianity of a state 
increases by increasing either the intensity of the input states or the conditioning value in the CPS scenario. This last condition 
requires photon-counting detectors endowed with a good linear response, such as those we used in our experiment. 

The use of a class of diagonal states in the photon number basis (the thermal ones), allows us to obtain a high degree of control 
between the analytical theoretical expectations and the experiment, which is a relevant point in view of further investigations. In 



particular, we are planning to apply our schemes to more exotic classical states, such as the phase-averaged coherent states M20I1 



These are characterized by a non-Gaussian nature themselves, and, thus, the possibility to perform conditional, non-Gaussian 
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FIG. 1 1 : (Color online) Reconstructed photon-number distributions for three different examples (black, green, red) of unconditional states 
(full symbols) and for the corresponding IPS states (empty symbols). The theoretical curves are plotted as lines according to the same choice 
of colors. The fidelity / is also reported. 
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FIG. 12: Log-linear plot of the lower bound s[^ips] for the non-Gaussianity <5[f?ips] as a function of the total mean detected photons Mt + Mr: 
experimental data (dots) and theoretical curve (line). 



measurements on them becomes particularly intriguing. Moreover, in this case analytical calculations may be carried out only to 
a certain extent: this is a clear example in which the reliability of the setup is a key point, as theoretical expectations are limited 
to numerical results. 

Though we only focused on classical states, our experimental procedure opens the way to further developments toward the 



generation and engineering of more sophisticated quantum states by considering non-classical states as the inputs 1 22, 



28D, which 



may be useful for quantum information protocols involving non-Gaussian states, such as entanglement distillation protocols 
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Appendix A: Experimental lower bound for the non-Gaussianity 



For a single-mode state diagonal in the Fock basis, i.e., g = J2 n Pn\n)(n\, the nonG measure [25] is given by 

S[g] = S[u(N)] - S(g) = S[u(N)] + logp„, (Al) 

n 

where v(N) is a thermal state with mean photon number TV = J2 n n Pn- Being based on the knowledge of the actual photon 
distribution p n , the calculation of S[g] requires measuring with an ideal (i.e., with unit quantum efficiency) photon-number 
resolving detector. In the presence of inefficient detection, one can only retrieve the detected photon number distribution q m = 
Tt[f>II m (??)], where H rn (rj) is given in Eq. (0 and rj is the quantum efficiency. Nevertheless, in the following we will show that 
the quantity 

e[g] = S[v{M)] + ^q m logg m , (A2) 

m 

where M — J2 m m 1 m ~ V-^, is a lower bound for the real non-Gaussianity S[g], i.e., e[g] < 6[g\. Note that since e[g] depends 
only on q m , it can be calculated starting from the experimental results. 

The inefficient photodetection process can be described by mixing the quantum state g with the vacuum at a BS having 
transmissivity r\ followed by perfect detection on the transmitted beam, thus obtaining 

q m = Tr[£( Q )\m)(m\], (A3) 

where £(g) = Tt2[Ubs(ii)q ® |0)(0|E/jg S (T7)] is the lossy Gaussian channel. Since g is diagonal in the Fock basis, £(g) is still 
diagonal 

£(e) = X>« s (l n K n l) = E^l m >H, (A4) 
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in which q m = X^m^™- 8 ™,™^)- To obtain £{g) we used £(\n)(n\) = Y%=o B n,l(v)\l){ l \> with B n,i(v) defined in Eq. ©. 
By using Eq. (I Alt , we obtain 



S[£(g)] = S[u(M)} + Qn \ogq n = e[g}. 

n 

As the nonG measure S[g] is non-increasing under Gaussian maps 12511 . we finally get 



(A5) 



e[g] = S[£(q)] < 5[g}. 



(A6) 



Summarizing, given a quantum state g, diagonal in the Fock basis, we can measure the probability distribution of the detected 
phtotons q m and evaluate Eq. JA2b as a lower bound for the actual non-Gaussianity 6[g\. 
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